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ABSTRACT 

The Teerikorpi incompleteness bias in the distance modulus of a galaxy cluster 
that is determined from incomplete data using the Tully-Fisher (TF) method 
is discussed differently than has been done in earlier papers of this series. A 
toy cluster is made with zero intrinsic TF dispersion but with slopes that differ 
between the calibrators and the cluster data, showing the bias caused by incorrect 
slopes. Intrinsic dispersion is added to the model and two strategies are used to 
analyze the data; first by binning the data by line width and then by apparent 
magnitude (the direct method), and second by binning by magnitude and then 
summing over all line widths (the inverse method). To illustrate these strategies, 
a composite cluster is made by combining the observations of Virgo A and B 
subclusters with those for the Ursa Major I and II clusters, corrected to the 
Virgo A distance. The cluster data are calibrated using Cepheid distances to 
25 galaxies that have adequate TF properties. Different moduli calculated with 
varying completeness limits are displayed. The cluster modulus derived from the 
complete cluster sample gives (m — M)° = 31.42 ± 0.2 (external) for Virgo A, 
31.80 ± 0.16 for Virgo B, 31.26 ± 0.13 for UMa I, and 31.58 ± 0.17 for UMa II. 
Combining the Virgo A distance (D = 19.2 Mpc that has a range from 17.5 to 
21.1 Mpc) with its expansion velocity of 1175 ± 50kms _1 as tied to the remote 
kinematic frame gives a Hubble constant of 61 km s _1 Mpc -1 with a range from 
53 to 70. 



Subject headings: Research — Tutorial — Galaxies — Data Analysis and Tech- 
niques 
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INTRODUCTION 



The literature on the effect of bias in the measurement of extragalactic distances is 
extensive. Taken is a whole it seems comprehensive and exhaustive. Why then propose 
again an expose of parts of the subject? The problem is not so simple in different situations, 
and new illustrations of the details can be useful. 



A significant review, complete to 1996, is by ITeerikorpil (119971) . He treats many aspects 
of the probl em beginning with the early insigh t s of iKaptevnl (119141 ). and then proceeds to 
the classical iMimauistl Jl92p| . Il922h and lScottl (|l957j) effects. 



The Malmquist calculation only concerns the bias integrated over a total sample. It 
does not show how the biased distance errors change as the truncation level is made fainter. 
Therefore, the classical Malmquist bias is not useful in correcting individual distances at 
various magnitude levels above the truncation level. For that, more complicated correction 
procedures mentioned below, are needed. 

An equally complicated problem is to analyze the bias effects for individual galaxy 
distances using the Tully-Fisher line width (LW) method applied to field galaxies chosen in 
different ways, such as from magnitude-limited catalogs, or by imposing apparent diameter 
limits on 21-cm lists before measuring the radio line widths, or by restricting line widths at 
either the high and/or low LW limit, etc. 

The literature on the use of the 21-cm line as a distanc e indicator beg a n with the dis- 
covery of th e correlation of LW w ith absolute magnitude by [Roberts! (119691 ) . iGouguenheim 
( 119691 ). and iBottinelli et al.l (119711) . The correlation was then developed as a distance indi- 
cator bv lTullv & Fisheri Jl977T ) 



The bias properties of individual TF distances for field galaxies, chosen from magnitude- 
limited, or apparent diameter-limited samples, is clearly summarized by Teerikorpi in his im- 
portant review. Cor rection procedures for b i as can be made either by (1) the method of "no r- 
malized distances" ITeerikorpil 1 19841 . Il99pl : IBottinelli et~aD Il986l : iTheureau et~aD Il997afl ) 
where a "plateau" is seen in an M, (m — M) diagram that is calculated using TF distances 
and when all galaxies not on the plateau are dis carded, or (2) corrected by a method such 



as the "triple ent ry" procedure 



Sandage 



Sandage. Tammann. fc Federspiel 



1994a 



hauer diagrams (Sandage! 1988a bT 1994a bT 1999 ; Federspiel. Sandage &: Tammann 1994 



b) or by the equivalent method of Spaen- 



1995 



The problem of bias for cluster galaxy samples using the TF method would seem at 
first to be simpler than for TF field galaxy samples. Occasional early comments in the 
literature stated that bias problems for clusters do not exist because all members of the 
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sample are nearly at the same distance. This is incorrect . Bias still e xists when a cluster 
luminosity function is sampled incompletely. iTeerikorpil (119871 . Il990l ) has called this the 
"cluster population incompleteness bias", of which there is now a significant literature. 



A listing of the central pap ers on the incompleteness bias include those by iBottinelli et al 



f l987 ). lBottinelli et a 

Jl990h. IWillick! Jl994f). IFST94I (their Figs. 5, 6, fc 8), 



( 1988a b), Kraan-Korteweg. Cameron, fc Tammann ( 19881 ). Fouque et al 

STF95I (theixFigs. 2, 3, & 10), 



Giovanelli et al.l (ll997al Jbl). iTeerikorpi et al.l (11999). Masters et al.l (20061) althoug h we dis 
agree with part of their analysis of the bias (ITammann. Sandage. fc Reindll 120071 ). and by 



others cited therein. Why then is another paper on the cluster incompleteness bias useful, 
and what is the justification of this one? 

Clarifications can be made on several points that include these. (1). Bias will exist with 
incomplete cluster sampling due to two separate effects, and it is useful to separate them. 

(a). Applying an incorrect TF slope from the calibrating galaxies to the cluster sample 
is the principal reason for bias. Only when the adopted slope for the calibrators is the same 
as that of the sample can a slope bias be eliminated even if the cluster luminosity function is 
sampled "completely" in a certain way. This is true both for the "direct" and the "inverse" 
TF formulation, i.e. reversing what is the independent and dependent variables in the least 
squares regressions of LW on apparent magnitude. 

The direct and inverse slopes will always be different in an y correlation of two indepen- 
dent variables if there is intrinsic dispersion in the correlation (cf. lSeareslll944j ; lFeigelson fc Babu 
19921 ). 



For the Tully-Fisher correlation it turns out to be crucial to use the correct slope in 
analyzing the data in either the direct or the inverse mode. However, there is the question of 
whether we should use the direct slope, the inverse slope, or something in between, sometimes 
called the compromise slope, in comparing the observed TF relation with any particular TF 
calibration. 

(b). The second but related problem is the effect of the intrinsic dispersion of the TF 
correlation itself in the presence of incomplete sampling, even if a "correct" slope is used. 
The in yerse formulation is free from bias in a sampling to d i fferent apparent magnitude 
limits (jSchechterl 11980c iTullyi Il988l ; iHendry fc Simmons! 1 1994 ; lEkholm fc Teerikorpil 119971 ; 
Teerikorpi et al.lll999l ) if the LW distribution is complete (not truncated). However, this is 
true only if the inverse slope is used for both the calibrators and the cluster sample rather 
than the more usual practice of using the direct (or the compromise) slope that has been 
determined from the sample but which is then imposed upon the calibrators (§ H] below) in 
using the inverse method for the sample. Rather, the inverse slope must be used for both 
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the calibrators and the sample. 



(2). It has sometimes been written that the incompleteness bias would disappear if 
the scatter in the observed TF cluster LW-magnitude correlation (i.e. either an intrinsic TF 
dispersion or a back-to-front depth effect, or both) would be zero. This would only be true 
if the correct slope for the calibrators is used. However, such a slope is never known because 
of errors in the calibrator distances, giving a dispersion in the TF calibration solutions. This 
is the reason for the difference between the direct and inverse least squares solutions for the 
calibration. One of the purposes of this paper is to explore the effect of this dispersion both 
in the calibration and in the galaxy sample using a variety of TF slopes. 

The effect of calibrator dispersion on the derived TF slope, especially in the direct mode, 
can be made particularly transparent using real calibration data from the HST Cepheid 
database for TF calibrating galaxies. We set out in § H] a new Tully-Fisher calibration using 
our HST Cepheid distances (ISaha et al.ll2006l ). These are the basis o f our 2006 distance scal e 
with its determination of the Hubble constant (ISandage et al.ll2006l ; iTammann et al.ll2007l ). 



(3). A different formulation of the TF bias problems can be made and c an be compared 
with earlier discussions that may seem to have been unnecessarily opaque (eg. ISandagdll994bl . 
Paper II). The present paper aims to clarify the TF bias problem, recasting the discussion 
in Paper II of this series by retracing the steps this author recently took in considering the 
subject again. The purposes here are fivefold. 

(A). A simplistic toy (zero intrinsic dispersion) model is set out in § [2] showing what 
happens to the calculated distance moduli when an incorrect slope is used for the Tully- 
Fisher correlation for either the calibrators or the cluster data. A restriction to avoid a bias 
is stressed requiring that the depth of the magnitude sampling must be faint enough to be 
symmetrical about a "crossing point" in the B^-UN correlation between the direct and the 
inverse regressions if the slope is incorrect. 



(B). A more complicated toy model with intrinsic dispersion is discussed in 
and without the correct slopes for the direct and inverse regressions. 



[3] with 



(C) . Real data are used in § H] to make a mock (composite) cluster by combining ob- 
servations for the spirals in the Virgo A and B cluster cores with similar data for the Ursa 
Major I and II spiral aggregates, corrected for the small distance differences. 

(D) . §0sets out a new calibration of the TF relation us ing our HST metallicity-corrected 
Cepheid distances (ISaha et al.l 120061 ; ISandage et al.ll2006l ). This calibration using various 
slopes is applied to the composite cluster for both the direct and inverse regressions in § [6j 



(E). The resulting distance to the Virgo Cluster A spiral core is set out in §[7J The Hub- 
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ble constant derived from this distance, using the cosmic expansion velocity of the cluster, 
freed from all local non-cosmological flows by tying to the distant cosmic kinematic frame, 
is set out there. 



2. MODEL OF THE BIAS CAUSED BY INCORRECT TF SLOPES FOR 
THE CASE OF ZERO INTRINSIC TF DISPERSION; THE 
IMPORTANCE OF SYMMETRICAL SAMPLING OF THE 
MAGNITUDE INTERVAL 

As said, even if the intrinsic TF dispersion (and/or an appreciable depth effect) would 
be zero, there would yet be a bias in a derived mean cluster distance if (1) an incorrect 
slope for the TF correlation is used, and (2) if the sampling is either incomplete at faint 
magnitudes or is non-symmetrical about the cross-over point of the direct and inverse slopes 
in the TF diagram. We illustrate these points using the simplistic model shown in Figured] 
that has a zero TF intrinsic dispersion. 

The insert in Figure [1] shows a dispersionless TF correlation for an imaginary cluster 
that has a true slope of dm/d log u^o = — 7 and that is normalized at B^ 1 = 11.5 for 
logu^o — logW-y sin i = 2.50. A calibration of the TF relation with a slope of —7 that 
gives (m — M) = 31.5 for the toy cluster is in the main body of the diagram, normalized at 
absolute magnitude —20 for logti^o = 2.5. Two deviant slopes for the calibration are shown 
with slopes of —9 and —5. These cross the true slope at M = —20 and logu^o — 2.5 by 
construction. These imitate (albeit with exaggeration) the inverse and direct least squares 
regressions for real calibration data (Figs. [5] and E] in § 14. ip . For real clusters the inverse 
correlation is always steeper than the direct. (The actual slopes for the 25 galaxy Cepheid 
calibration in §|6]are —6.161 for the direct regression and —8.475 for the inverse, from Table[7] 
later) . 

We can illustrate the bias on the calculated cluster moduli by using each of the cali- 
bration (incorrect) slopes of —9 and —5. Consider first the calibration slope of —9 applied 
to the toy cluster. This steep slope gives too bright an estimate of the absolute magnitude 
for all line widths larger than the cross-over point at logu>20 — 2.5, and hence the distance 
moduli are too large for all logti^o > 2.5 and too small for line widths that are smaller than 
logu^o = 2.5. The opposite is true for a slope of —5; galaxies with logw 2 o larger than 2.5 
will be incorrectly calibrated too faint, giving smaller (incorrect) distances for logu;2o > 2.5 
and too bright giving larger (incorrect) distances for logu^o < 2.5. 

In a complete sampling for LWs above and below logu>2o = 2.5, and if the number of 
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galaxies along the line in the insert diagram is constant, then the average distance modulus for 
the complete sample will be correct at (m — M) = 31.5. But this is true only if the summing 
over all apparent magnitudes is carried out for magnitude limits that are symmetrical about 
the cross-over point at logw^o = 2.5 (or B^ 1 = 11.5 in the toy model). Said differently, a 
correct mean modulus will only be obtained by summing over magnitudes that go to equal 
limits above and below 11.5, i.e. say from 10.0 to 13.0. If the magnitude boundaries are 
not symmetrical about this cross-over point (here at Bji 1 = H-5), or if there is a population 
gradient with magnitude, then the derived mean modulus will be incorrect. Of course, also 
if the summing over magnitude is incomplete (say only over the restricted interval of 10 to 
12 mag), an incorrect mean modulus will be derived if the wrong slope is used. 

This is the classical "population incompleteness" bias, but we see here that it is more 
complicated than simple incompleteness. It depends on slope errors and the necessity for 
symmetrical sampling above and below the cross-over point of the direct and inverse least 
squares solutions even if there is no gradient in galaxies numbers along the ridge line. This 
crucial point (see §Ej) is illustrated in Tabled] for the dispersionless toy cluster that has been 
put at a true distance modulus of (m — M) = 31.5. 

The first half of the table gives the calculated modulus at the fixed magnitudes in 
column (1) using arbitrary TF slopes that range from —9 to —5. Fixing all cross-over points 
of these lines to be M = —20 at logw^o = 2.5 gives the equations of the assumed lines in 
the main body of the diagram as M = — 91ogu> 2 o + 2.5 to M = — 51ogw 20 — 7.5, etc., and 
B^ 1 = — 71ogw 2 o + 29.00. All entries for (m — M) at the stated apparent magnitudes in 
column (1) of Table [T] follow in an obvious way from these equations. 

The calculated moduli for the slopes of —9 and —8 with apparent magnitudes brighter 
than 11.5 (the cross-over point) are larger than the true modulus of 31.5 and visa-versa for 
fainter magnitudes. The opposite trend exists for slopes of —6 and —5. These results are 
shown in Figure[2] (top). Of course, for the "true" slope of —7, all calculated moduli are 31.5 
by construction. 

Using the calculated moduli that are valid at the specified apparent magnitudes, we 
can sum the entries to obtain the average modulus value for all galaxies that are brighter 
than any particular limiting magnitude into the cluster luminosity function, shown in the 
bottom panel of Figure [2j These moduli are, of course, those that would be obtained in a 
real case by averaging the modulus values for all galaxies in the sample that is complete only 
to brighter magnitude limits. This is the incompleteness bias. 

The effect of an incorrect slope on the bias using incomplete sampling is set out explicitly 
in the second half of Table [U and shown in the bottom panel of Figure [2j The points to 
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note are these. 

(A) . For steeper slopes than the "true" slope (i.e. for —9 and —8), the derived mean 
distance moduli decrease as the sampling is done deeper into the cluster luminosity function. 
The bias effect from B = 10 to 13 is 0.4 mag in the modulus differences for the slope of —9 
and 0.2 mag for the slope of —8. 

(B) . The opposite is true for the shallow slopes of —6 and —5. 

(C) . The correct modulus is obtained only by sampling to B = 13. This is three 
magnitudes into the cluster luminosity function. 

(D) . Sampling fainter than B = 13 gives incorrect modulus values (except for the —7 
case using the "true" slope). This is the effect of "non-symmetrical" sampling about the 
cross-over point. 

Of course, this toy model is too simplistic because there is no dispersion in the TF 
calibration nor in the TF correlation for the cluster galaxies. 

3. A MODEL OF THE BIAS DUE TO AN INTRINSIC TF DISPERSION 
AND/OR INCORRECT SLOPES FOR THE TF RELATION 

To understand the bias properties of the cluster TF data it is useful to analyze the 
LW-apparent magnitude data in two ways. The first is to bin the data into intervals of u>20 
and then, within each interval, to order the listings by apparent magnitude. The second is to 
bin the data by apparent magnitude and then to order by line width within each magnitude 
interval. 

The first way is illustrated in Figure [3] where the intrinsic dispersion (plus any back-to- 
front variation) is shown by the dotted envelope lines placed symmetrically about the ridge 
line and where three line-width intervals are drawn for illustration. An arbitrary magnitude 
limit is put at m = 11.5, which is about 1.5 mag fainter than the brightest galaxy in the 
cluster in this example. 

All galaxies in the line-width interval farthest to the right are brighter than the magni- 
tude cut-off, whereas galaxies at smaller LWs to the left are progressively lost to the sample 
at smaller and smaller LWs. This causes the bias. 

Consider first the unbiased LW interval to the right. Due to the dispersion, some 
galaxies are brighter than the central ridge line and some are fainter. If the position of the 
ridge line has been calculated by a least squares regression using magnitude residuals at fixed 
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W20, i-e. in the "direct" solution, there will be equal numbers of galaxies above and below 
the calculated ridge line. This line is the most probable apparent magnitude at that w 2 o, 
although, because of the dispersion, no galaxy will have the ridge-line apparent magnitude 
unless, of course, it is on the ridge line. 

Consider next how the data in this LW-interval are used to obtain a system of distance 
moduli using some absolute TF calibration such in Figure [8] later in § [51 The calibration 
will, itself, have a dispersion that will be a combination of the intrinsic dispersion due to 
the physics that governs the TF correlation and the measuring errors in the distances of the 
calibrating galaxies. 

When the ridge-line calibration (i.e. the central line in Fig. [8] later) is applied to the 
ridge line of the observed cluster TF relation such as in Figure El a modulus value is obtained 
of each galaxy in that w 2 o interval. But, as said before, the calculated modulus of any given 
galaxy so obtained is incorrect, unless it is on the ridge line. 

In detail, consider the unbiased logu^o interval from 2.70 to 2.75 for the right-hand 
strip in Figure El Suppose that all galaxies in the cluster are at the same distance (no 
back-to-front effect) and there are no errors of measurement either in W20 or magnitude. 
For this case, galaxies near the upper envelope line are those that are actually overluminous 
for their line width due to an intrinsic dispersion (i.e. from the physics) of the TF relation. 
Hence, applying the most probable calibration (which is fainter here than these particular 
galaxies) from the ridge-line calibrator relation to such intrinsically bright galaxies will give 
too small a calculated distance modulus for them, and visa- versa for galaxies that are actually 
underluminous (i.e. below the ridge line). Nevertheless, if the strip is filled symmetrically 
above and below the ridge line, and if the slope of that line has been determined by the 
"direct" least squares regression, then the mean modulus found by averaging the individual 
moduli in the strip (all of which are incorrect except for those on the strip) will be the correct 
modulus, to within statistics. 

However, if the strip is not sampled completely, as in the second and third strips in 
Figure [3J where a magnitude cut truncates the distribution, the average of the individual 
moduli will be progressively in error at the smaller line widths as LW intervals move toward 
the left. The average of modulus values for galaxies above the magnitude limit line will then 
be too small by amounts that will decrease as the fraction of the accessible sample (the ratio 
to the complete sample) increases as the magnitude limit lines are moved fainter. This is 
the incompleteness bias. 

Nevertheless, averaging over all line widths and to a sufficiently faint apparent magni- 
tude will give the correct distance modulus if, of course, the strips are filled symmetrically 
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above and below the most probable ridge line. 

This bias effect produces an error in the calculated slope determined from incomplete 
data in the direct regression. The error will be a function of the depth of penetration into 
the cluster luminosity function. This is obvious from Figure |3j The slope appropriate for the 
unbiased LW interval at the right (logw^o between 2.70 and 2.75) is the correct direct slope 
(except for statistics) that would be determined from the complete sample using the complete 
luminosity function. However, if we use only data brighter than any brighter magnitude cut, 
here put at 11.5 mag, the midpoint average magnitudes for the biased LW intervals (all to 
the left of the unbiased interval in Fig. [3]) are all brighter than the true ridge line (unbiased) 
points shown by the three dots in Figure [31 Hence, the slopes that would be determined 
by calculation using only data brighter than any magnitude cut-off, will all be too shallow, 
approaching the correct slope only when the magnitude limit is so faint so as to include the 
complete sample. 

Hence, because of the incompleteness bias, the calculated mean modulus using the direct 
slope should vary pro gressivel y with cut-off magnitude. This is the explanation of why in 



our earlier discussion (ISTF 95l . Figs. 2, 3, 6, and 10) the calculated slope of the observed TF 
correlation in the direct formulation varies with cut-off magnitude. This effect on the slope 
will be demonstrated in § 14.11 (Fig. [7]) using actual data. 

The ideal schematic model of Figure [3] is still too simplistic because the effect of using 
an incorrect slope for the ridge line has been ignored. We mean by this that the (direct) 
slope shown in Figure [3j as it would have been derived from the complete sample, is assumed 
to be without error. 

With this in mind, consider the inverse correlation where the calculated TF slope is in- 
dependent of the depth sampled into the cluster luminosity function, in contrast to the direct 
slope calculated from incomplete data. The strategy of analysis for the inverse calculation 
is shown in Figure HI 

In Figure IU as in Figure [3J the scatter of the cluster data at a given apparent magnitude 
is shown for a TF correlation that has intrinsic and/or back-to-front scatter. The inverse 
ridge line (marked I) is made to go through the midpoint of each shaded magnitude interval, 
because that is what the least squares regression gives by its procedure of using residuals in 
the LWs at given magnitudes rather than visa versa. 

If the distribution of W20 within each magnitude interval is symmetrical about the ridge 
line, the average distance modulus of all the galaxies within each magnitude interval will be 
the true modulus to within statistics. This average at every magnitude interval will be bias 
free. However, as is evident in Figure HI this will only be true if the inverse slope is used. It 
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will not be true if the direct slope is used with the inverse method. 

Using either the I or the D ridge lines, and imposing in turn their slopes on the calibrator 
sample, gives individual modulus values for each galaxy in any of the individual strips in 
Figure HI Again, as before, these individually are incorrect, except for galaxies precisely on 
the ridge line. Nevertheless, if the distribution of modulus values within each magnitude 
interval is symmetrical about the "I" ridge line, then the average of the individual modulus 
values within each magnitude interval will be the true value provided that the I slope has 
been forced on the calibrators themselves. Clearly, as seen from Figure HJ the same is not true 
using the direct slope. There are more galaxies to the small LW side of the TF distribution 
for magnitudes brighter than the cross-over point, and fewer galaxies to the left of the D 
line for fainter magnitudes than the crossover point. This produces the bias in the average 
modulus using the direct slope for all magnitude cut-offs brighter than the limit at B — 14. 
This is the cluster incompleteness bias in the direct formulation if the sampling is incomplete. 

In the inverse method we need only average all the modulus data in a given magnitude 
interval to get the correct cluster modulus, in contrast to the direct method where we need 
to average over the whole of luminosity function, or we must apply the correction methods in 
the direc t formulation de rived either from the "normalized distance" method of the F rench 
workers (ITeerikorpil 119971 ). or by the method of Spaenhauer diagrams JSandagelll994al lbll 



Figure H] also shows the effect caused by an error in the adopted slope in the inverse 
method. This will cause a bias in the derived distance for the reason just described of using 
the direct slope when the inverse slope should be used. 

Consider this effect of using the direct slope with the inverse method and summing the 
data only to a given apparent magnitude. Consider first all magnitude intervals that are 
brighter than the cross-over point at m — 11.5 and logu^o = 2.5. An average modulus 
that is too small will be calculated at each of the bright magnitude intervals because there 
will be more galaxies to the left of the adopted "D" ridge line than to the right (i.e. giving 
moduli values that are too small compared with the ridge line average). The opposite will 
be true for all magnitude intervals that are fainter than m = 11.5 (seen also in the top part 
of Table [T]for a toy cluster with no dispersion). Only by summing the total data over all 
magnitude intervals can the correct true modulus be obtained. However, as in Figure [1] and 
Table [lj the sum over all magnitudes in a complete sample must be made no fainter than 
the magnitude level that is symmetrical about the cross-over point, showing again the need 
to adopt a symmetrical magnitude limit about the cross-over point in both the direct and 
inverse solutions. 



Hence, on both accounts of the effect of (1) incomplete sampling using the direct TF 
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slope as in Figure [3] and then summing over all magnitudes (usually such completeness is 
not available), or (2) of summing over LW to given magnitude limits (Fig. H] for the inverse 
method), the conclusion is that the true inverse slope must be used to avoid both (a) the 
error in the adopted slope and (b) the incomplete sampling effects. These conclusions will 
be demonstrated in § [6] using real data for a composite (artificial) cluster. 



4. A COMPOSITE CLUSTER MADE TO ILLUSTRATE THE 

INCOMPLETENESS BIAS 



4.1. Construction of the Composite Catalog 

To demonstrate the points made in the last two sections it is useful to analyze real LW 
and magnitude data for a galaxy cluster. Only a few catalogs exist that are complete enough 
for real clusters to be useful for such a demonstration. The most complete are for the Virgo A 
and B sub clusters and for the U rsa Major "clust e r", ea rly called the Ursa Major Cloud by 
Hubble and identified as such by lHumason et ajj (119561 ) in the Mount Wilson redshift list. 



However, the number of spirals that are suitable for a TF analysis (inclinations greater 
than 30° and adequate photometry to derive fully corrected apparent magnitudes) in these 
three aggregates is so small that small-number statistics can mask the bias effects we are 
seeking. To increase the size of the sample we have created a composite cluster by combining 
the data for the Virgo A and B subclusters with those for the Ursa Major cloud to form a 
single catalog. Corrections have been made to the apparent magnitudes of the galaxies in 
the Virgo B cluster and the Ursa Major Cloud to compensate for the small differences in the 
distances of each aggregate relative to Virgo A. 



The data for the Virgo subclusters are taken from the listings by lFederspiel. Tammann. fc Sandage 



( 119981 ) in thei r Table 3. Membership for Virgo A and B are based on the Virgo Clus- 
ter Catalog of iBinggeli. Sandage. fc Tammann! (119851 ) as reanalyzed for cluster members by 
Binggeli. Popescu. &: Tammann! (Il993l ). Only galaxies with inclination greater than 30°, and 
only those assigned either to subclusters A and B are used. There are 43 galaxies in Virgo A 
and 16 in Virgo B in the composite catalog. 



The Ursa Major Cloud has been studied in the p ast by ISersid (Il960f). Ide Vaucouleurs 



1975 . especially his Fig s. 8 & 9) in his review of groups, iTullv fc Fisherl Jl987h 



Verheiienl ( 1997 , his Fig. 7 in chapter 4 of his thesis) , Tully et al. (jl996 ) , 



'ierce &: T ullv 



Peletier fc Willner 



Federspiell (119991 ). and undoubtedly others. 



We use here the cluster membership of iFederspiell from his Table 3.2 for the spirals that 
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are suitable for the TF method where he lists the necessary B^ 1 , line width, and V220 radial 
velocity d ata (t o test members hip). Federspiel's membership criteria, similar to those of 
Verheijenl and of lPierce fc Tulfyl . are (a) heliocentric redshift interval from 700 to 1500 km s _1 , 
and (b) angular distance within 7.5° of RA(1950)= ll h 54 m , Dec(1950)= +49.5°. 

The radial v elocity distribution in 1*220 , corrected for Virgo infall by the model of 
Kraan-Kortewegl ( 119861 ) . shows a bimodal character with means at 1060 ± 29kms _1 and 
1485 ± SSkms™ 1 , indicating an appreciable depth effect along the fila ment connect i ng th e 
cloud with the Virgo complex that is shown so remarkably in Figure 6 of lKlypin et al.l (120031 ) . 
We treat separately the two distributions, assigning galaxies with V220 < 1250 km s -1 to be 
members of an Ursa Major I aggregate, and those with V220 > 1250 km s -1 as members of 
Ursa Major II. There are 18 spirals suitable for the TF method in UMa I and 15 in UMa II. 

Hence, there are 92 galaxies in the artificial composite cluster (43, 16, 18, and 15) from 
Virgo A, Virgo B, and UMa I and UMa II. 

The B^ 1 magnitudes used here, corrected to a "total" magn itude, are within a few hun- 
dredths of a magnitude of the RC3 (jde Vaucouleurs et al.lll99ll ) system. (Slight differences 
at the few hundredths mag level exist because no K redshift corrections were applied here, 
differing from the RC3, and some updating was mad e) . The L W data for the four aggregates 
are on the homogeneous line-width system used by IFTS 98l . taken from the Lyon-Meudon 
Extragalactic Data Base (the LEDA). These LWs average 0.033 ± 0.002 dex smaller than 
the LWs that are calculated from the RC3 with the RC3 inclinations calculated from the 
log R25 values listed there. 

Small magnitude corrections between the four groups are needed to reduce all data to 
the distance of Virgo A. These were determined by comparing the zero points of the least 
squares TF correlations of Virgo A with those of Virgo B and UMa I & II. To calculate these 
corrections, the Virgo A TF slope is imposed on Virgo B and UMa I & II. These gave zero 
point ma gnitude differences such that the B^ 1 magnitudes for Virgo B are made brighter 



than the IFTS 98l listings by 0.38 ± 0.16 mag (Virgo B is more distant). The correction for 
UMa I (t he smaller velocit y aggregate) is 0.16 ±0.13 mag made faint er than the magnitudes 
listed by iFederspiell (Il999l ) (UMa I is closer than Virgo A), and the iFederspiell magnitudes 
for UMa II (the larger velocity aggregate) were made brighter by 0.16 ±0.17 mag (UMa II is 
more distant than Virgo A) so as to put the magnitudes on the distance system as Virgo A. 

Using these magnitude corrections we have made an artificial cluster catalogue from 
which the TF correlation of log W20 with the B^ 1 magnitudes can be made. The result is in 
Figure [5] where all 92 galaxies of the artificial composite cluster are shown. A more select 
subsample of 84 galaxies is shown in Figure [6] where a magnitude cut is made at B^ 1 = 14.0, 
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eliminating the eight faintest galaxies. 

The slopes of the direct and the inverse least squares regressions with the cluster data 
are shown in parentheses in the diagrams. The equations of the regression lines are: 

£°'*(RC3) = -6.622 log w 20 (LEDA) + 28.195 ± 0.103, (1) 

for the direct regression using all 92 galaxies in Figure \5[ and 

5^(RC3) = -8.032 log w 20 (LEDA) + 31.598 ± 0.114 (2) 

for the inverse correlation for the same sample. 

The equations for the more restricted 84 galaxy sample (cut at B^ % = 14.0) are: 

B$\BC3) = -5.936 log w 20 (LEDA) + 26.453, (3) 

for the direct regression, and 

£^(RC3) = -7.547 log w 20 (LEDA) + 30.391. (4) 

4.2. The Sensitivity of the TF Slope with Sampling Depth into the Cluster 

Luminosity Function 

The data for the artificial cluster have been binned by apparent magnitude, enabling 
a study of the effect on the biased slope of the Tully-Fisher correlation due to incomplete 
sampling. Least squares solutions have been made from the composite cluster data as the 
magnitude cut is put progressively fainter. The results are listed in Table [2] in columns (2-6) 
for the total sample of 92 galaxies to the limit of B = 15, and for the restricted sample of 
78 galaxies in columns (7-9) to be discussed later in § [61 

Column (1) shows the magnitude cut into the luminosity function of the composite 
cluster. The number of galaxies to each truncated limit is in column (2). The slope and 
zero point for the direct regression (residuals taken in magnitude at a fixed LW as in m = 
alogw+b), are in columns (3) and (4). The slope and zero points of the inverse correlation, as 
in logw 2 Q = cB^ 1 + d, are in columns (5) and (6). The slopes and zero points for the complete 
sample of n = 92 for B > 15 (the last row) are, of course, the same as in equations (pQ) and 
d2J). Columns (8) and (9) show the calculated least squares direct and inverse slopes for the 
restricted sample (n = 78) to the listed cut-off magnitudes. The number of galaxies to that 
magnitude limit is in column (7). 



-14- 



Figure [7] shows the data in Table [2] where the progressi ve increase of slope for the direct 



correlation, consistent with Figs. 2, 3, 6, and 10 of ISTF 951 . is evident. The top panel shows 
the increase of slope with increasing faintness for the direct formulation using the complete 
sample of 92 galaxies (col. 3 of Table [2]). Less clear but consistent within statistics, is the 
lack of a correlation of slope with faintness for the inverse correlation. The bottom panel 
shows the same effect more clearly for the restricted sample of 78 galaxies from columns (8) 
and (9) of Table [3 

The expectation from § [3] (eg. from Fig. [3]) is that there should be a progressive change 
of slope as the cluster is sampled more deeply in magnitude using the DIRECT regression 
(residuals in magnitude at given LW), whereas no systematic variation should be present 
in the INVERSE solution (residuals in LW at given magnitudes). Figure [7] shows the slope 
effect in the DIRECT (dots) formulation. It is consistent with Figures 2, 3, 6, and 10 of 



STF95 



The situation is less clear in support of the prediction of no-slope variation for the 
inverse regression (open circles) because of the dip between 12 and 13 mag and the increase 
at fainter magnitudes. However, inspection of the distribution of galaxies within the envelope 
lines in Figures [5] and [6] shows the reason is the statistical noise caused by fluctuations in the 
small number statistics of the actual data, and the evident non-uniform filling of the vertical 
and horizontal bins i n Figur es [3] and HI seen by inspection of Figures [5] and [6] here, and 



especially Figure 7 of IFTS 98l . In Figures \5\ and O for the magnitude interval from B = 12 to 
13 there are more galaxies to the left of the ridge lines than to the right, and conversely, from 
magnitude 13 to 14 the opposite is true. This causes the entries in the later Tables [5] and [61 
used to construct a magnitude- modulus diagram similar to Figure [2], to show this statistical 
noise as curvature from B = 12 to 14, reducing the ideal symmetry of Figure [2] for the toy 
cluster to the noisy reality of the real data for the composite cluster. We shall encounter 
the problem again in § [HI where the data for the composite cluster using various TF slopes 
resembles Figure |2j This statistical noise between B — 12 to 14 degrades the predictions of 
the ideal noiseless toy cluster of Figures [1] and [2j 

To this point we have not needed an absolute magnitude calibration of Figures and [H] 
to study the bias problem. In particular we have not needed a calibration to obtain Table [2] 
and Figures [71 but to carry the argument further we now do. 



- 15 - 



5. A NEW TF CALIBRATION USING HST CEPHEID DISTANCES FOR 

25 GALAXIES 

Equations ([UH]) and Figures [5] and [6] are only useful for determining distance moduli if 
a calibration of the TF regression in absolute magnitude is available from independent data. 
We use here a new calibration based on 25 galaxies in which Cepheids have been measured, 
the majority of have been observed with the HST space telescope. The new data are set out 
in Table [3j 

Column (1) lists the galaxy name. The type and luminosity class are in column (2) from 
the listings in A Revised Shapley Ames Catalog of Bright Galaxies (ISandage fc Tammann 



1981. 19871 ). Column (3) is the adopted distance modulus from ISaha et al. (bood . Ta- 



ble 8) based on period-dependent metallicity corrections requi red by the non-unique Cepheid 
period-luminosity relations that vary from galaxy-to-galaxy (Tammann. Sandage. fe Reindl 



2003 



Sandage. Tammann. fc Reindll 12004 ; iKanbur fc Ngeow 



2004 



Ngeow fc Kanbur 



2005 



Ngeow et al. 20051 ). These mod ulus values define our 2006 distance s cale based on type la 



supernovae and HST Cepheids (ISaha et al.ll2006l ; ISandage et al.ll2006l ). (The listed modulus 
for NGC4258 corrects the value listed in Table 8 of Saha which was a computational error). 



Column (4) is the fully corrected apparent magnitude on the scale of IFTS 98l which is 
close to that in the RC3. Column (5) is the M B 'Lx absolute magnitude which is column (4) 
minus column (3). The log line width at the 20% level and corrected for inclinat ion is in 
column (6) on the Lyo n-Meud on Extragalactic Data Base (LED A) as set out in IFTS 98 . 
Where no data exist in IFTS 98l we use the RC3 LW data made smaller by 0.033 dex. 



The difference in the distance moduli in column (3) with 16 galaxies in common with 



FTS 98l is 0.05 ± 0.03 mag with the distances in Table [3] being larger. 



Comparisons of the distance scale here with that of iFreedman et al.l (120011 . Cols 2 and 
5 of their Table 3) are these. The difference between Table [3] here and Table 3, column (2) 
of IFreedman et al.l (their old scale) is 0.09 ± 0. 04 mag from 25 o verlaps. Our distances are 
larger. The difference using the new scale of IFreedman et al.l (their Table 3, col. 5) is 
0.23 ± 0.04 mag, again our distances are larger. 

Comparing the larger sample of all Cepheids (not just those that make up the TF 
calibration) that are in common between our scale and Freedman shows (m — M) = 0.11 ± 
0.03 mag for Freedman's old scale (col. 2 of their Table 3) from 30 overlaps. Their new scale 
(col. 5 of their Table 3) compared with ours shows (m — M) = 0.25 ± 0.04 mag from 30 
galaxies in common. Again our scale is longer. 



The reason for the difference is that we account for the different slopes and zero points 
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of the Cepheid P-L rela tions from galaxy -to-galaxy according to metallicity corrections that 

are period dependent dSaha et al. 2006 ). The P-L slopes dif fer depending on metallicity 

(Fig. 5 and Table 4 of iTammann et al.1 120071 ). whereas those of iFreedman et all do not. For 
their new scale (column 5 of their Table 3), they adopt the shallow LMC P -L slope, not the 



steep er slope for the Galactic Cepheids which we have argued elsewhere ( ITammann et al. 



20071 ) that they should have used. 



The data from Table [3] are plotted in Figure M where the direct and inverse least squares 
regressions are shown with their slopes in parentheses. These slopes differ from the regres- 
sions for the composite cluster in Figures [5] and [61 They must not be used in analyzing 
the cluster data. Rather, the cluster slope must be imposed on the calibrator data (§ [6]) to 
obtain the calibrator zero point that is appropriate for the cluster data. This requirement 
is fundamental to avoid bias errors due to slope differences. This crucial point concerning 
the correct slope to use has sometimes been overlooked in the earlier literature on the TF 
method. 



6. BIAS IN THE DERIVED DISTANCE MODULUS DUE TO 
INCOMPLETE SAMPLING IN THE ARTIFICIAL CLUSTER USING 
BOTH THE DIRECT AND THE INVERSE TULLY-FISHER 

REGRESSIONS 

Individual distance moduli have been calculated for all galaxies (n = 92) in the artificial 
cluster using both the direct and inverse slopes. The calibration using the Cepheid data 
in Table [3] has been zero pointed by imposing the cluster slopes of —6.622 and —8.032 
from equations ([T]) and d2J) on the Cepheid calibrators. These give the Cepheid ridge-line 
calibrations of 

M°'j T) (FTS) = -6.622 log w 20 (LEDA) - 3.369 ± 0.103 (5) 
for the direct regression, and 

M°'j T) (FTS) = -8.032 log w 20 (LEDA) + 0.239 ± 0.114 (6) 

for the inverse. 

The individual (most probable, but all incorrect because of the dispersion except for 
galaxies on the ridge lines) moduli follow in an obvious way by combining these Cepheid 
ridge-line calibrations with the individual apparent magnitude and LW data for each galaxy 
in the cluster. 
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Following Figure [31 the data were first separated into log u>2o intervals in steps of 0.1 dex. 
They were then binned in apparent magnitude intervals of 0.5 mag and averaged. 

The average modulus in each W20 interval and each magnitude bin are shown in Table HI 
The number of galaxies that make up each average is in parentheses. 

The bias properties discussed in §[3] can now be illustrated using Table H] rather than by 
relying on hypothetical models. These properties and their consequences can most easily be 
understood by using Table @] and referring to Figures [3] and @] for the two binning strategies. 
There are three principal points. 

(1) . Table H] shows that within any given LW interval (columns 2-8) the average moduli 
show a systematic variation as the apparent cut-off magnitude is made fainter. This is 
expected from the discussion of Figure [3] in § [3J The variation is due to this. 

All galaxies in Figure [3] that are brighter than the ridge line magnitude at any LW have 
calculated moduli that are smaller than the ridge-line (the true) value and visa versa for 
those below the ridge line. Figures [Hk and [Hb show the trend of the individual moduli for 
both the direct and inverse slopes in the six LW intervals that range from 2.6-2.7 to 2.1-2.2. 

The expectation (§ [3]) that becomes manifest in Table HJA. is that the average over all 
magnitudes (i.e. summed vertically in each column of Table HJA.) should not vary system- 
atically with LW to within statistics using the DIRECT slope because the DIRECT ridge 
line threads strictly down the middle of each vertical column of Figure [3] by its least squares 
construction. This expectation is realized to within statistics by the averages set out at the 
bottom of Table H]A. for each LW interval. 

In contrast, using the INVERSE slope, the averages over all magnitudes at given LW, as 
analyzed via Figure HJ is expected to vary systematically with LW, meaning that the averages 
at the bottom of Table H)3 should show a strong systematic variation with LW because the 
INVERSE ridge line, not shown in Figure [31 does not thread the middle of the vertical LW 
columns. This expectation is well seen in the strong variation, outside the statistical noise, 
of the average moduli shown in Figure [Hb marked by the arrows, and seen in the bottom 
row of Table HE. 

(2) . However, rather than first binning by LW and then magnitude, the more realistic 
analysis, and the one always made in practice, is to work with averages made by summing 
over all line widths at various magnitude limits. The strategy is illustrated in Figure H] 
using progressively fainter magnitude cut-offs. Here, we are interested in the bias properties 
of Tables H]A. and B by summing over all line widths to various cut-off magnitude limits. 
The method using Tables H]A. and B is to sum horizontally (over all LW) in each magnitude 
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interval, and than sum over all magnitude intervals up to a given magnitude cut-off. 

The results are given in Table which is divided using the direct slope of —6.622 in 
the top half and the inverse slope of —8.032 in the bottom half. Column (1) shows the 
magnitude intervals used in Table [4] Column (2) is the number of galaxies used in the 
averages. Column (3) gives the average modulus for all galaxies summed over all LWs in 
this magnitude interval. These are calculated by summing the Table H] entries horizontally in 
each magnitude interval. Column (4) is the rms variation of these averages. Columns (5-8) 
give the average moduli calculated by summing the data in column (2) up to the magnitude 
cut-off listed in column (5). The number of galaxies making up this sum is in column (6). 
This is the accumulated number of galaxies by summing column (2). The mean modulus of 
the accumulated sums is in column (7). The rms of these averages is in column (8). 

(3). The expectation from §[3] is that the inverse slope must be used in summing hor- 
izontally in Table H] as in Figure H] to avoid bias. If so, there should be no systematic bias 
increase in the derived distance modulus based on the inverse slope as the grasp into the 
cluster luminosity function is deepened. However, this should not be true in the direct for- 
mulation where the derived average modulus is expected to increase with fainter magnitude 
cut-offs. 

These expectations are verified in Figure [10] where the data from Table [5] are plotted. 
Open circles are for the magnitude intervals in column (1). Crosses are for the accumulated 
data listed in column (7). The expectations from the discussion in §[3] are fulfilled to within 
statistics in the two panels of Figure [10J 

6.1. Effect of Using an Incorrect Slope for the TF Correlation 

Table [5] and Figure [10] are based on the formal least squares slopes for the TF data for 
the composite cluster for both the direct and inverse calculations. It is of interest to gauge 
the effect of systematic errors in these slopes by arbitrarily varying them over a wider range 
than given by the regressions on the actual. 

We have chosen TF slopes of —9, —7, —6, and —5 to add to the results in Tabled These 
cover the same range of slopes as in the toy cluster of Figures [T] and [2J The calculation was 
made by first redoing the Cepheid zero point calibration by imposing each of the assumed 
slopes on the Cepheid data in Table [3] These M, log u^o zero point calibrations were then 
applied to each of the 92 galaxies giving individual (m — M) modulus values for each cluster 
galaxy. These were then binned by apparent magnitude and summed over all LW for the 
galaxies brighter than a given apparent magnitude to produce a mean modulus for that 



-19- 



TF slope and apparent magnitude cut-off. The results are listed in Table [61 which can be 
combined with those in Table [5] for the slopes of —6.622 and —8.032 and Table [8] later for 
slopes of —6.346 and —7.519. 

Column (1) is the magnitude depth that is sampled into the cluster luminosity function. 
The number of galaxies making up the averages is in column (2). The average moduli to the 
cut-off magnitude are in columns (3), (5), (7), and (9). The rms of these listed moduli are 
in columns (4), (6), (8), and (10). 

Table [6] for the slopes of —9 to —5 combined with Table [5] for the data with the slopes 
of —6.622 and —8.032 show many of the general features of Figure [2] and Table [TJ 

As expected, the systematic run of moduli with apparent magnitude exists for all slopes 
that differ from the inverse slope of —8.032. For slopes steeper than this (—9 in Table [6]) 
the sign of the bias is that smaller moduli are calculated as the magnitude grasp into the 
cluster is increased. The opposite sign of the bias is evident (larger moduli for deeper cluster 
penetration) is seen for slopes flatter than —8.032. 

However, the detailed results, seen by plotting the data in Tables [5] and [6] (not shown), 
are not as clear as in the noiseless model of Table [1] and Figure [2j From Figure [2] (bottom) 
we expect that the null magnitude (where all modulus values converge) should be at B = 13, 
which is about 3 mag fainter than the brightest galaxy. However, the results for the composite 
cluster using real data show that the null magnitude is near B = 12 from Tableland between 
B = 12 and 12.5 from Table [51 As explained in § 14.21 the small-number statistical noise in 
the distribution of points between B = 12 and 14 is the probable cause of this difference. 
This shows the limitation on the accuracy of the TF method due to random noisiness within 
the dispersion when only a small number of cluster galaxies are available. Of course, the 
situation is improved using / magnitudes where the dispersion is smaller but the slope of 
the TF regression is near 10, adding additional uncertainties. 



7. DISTANCE TO THE VIRGO A SPIRAL CLUSTER CORE LEADING 
TO THE GLOBAL VALUE OF THE HUBBLE CONSTANT 

From the results concerning bias set out in the previous sections, how then do we use 
them to find the most probable distance modulus of the composite cluster? 

From analysis of the properties of the bias given above and from the results of Figures [TJ 
ITOl and [11] it might seem best to use only the moduli based on the inverse analysis because, 
given large enough samples, there is no bias. However, there are two reasons we have not 
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done this. 

(1) . The rms variation of the modulus summed to various magnitude limits, listed in 
Table and later in Tables [7] and El iri § 17.11 is always smaller for the direct formulation 
than for the inverse in a ratio of about 1.2 (eg. 0.576 to 0.488 from the penultimate rows of 
Table [5]). Hence, the direct moduli have greater weight than the inverse by a factor of about 
1.4. 

(2) . There are always differences in the final average moduli between the direct and 
inverse calculations that are independent of the bias problems. The cause is different Cepheid 
calibrations between the direct or inverse calculations. This is because a given galaxy will 
be assigned different absolute magnitudes depending on the two different calibrations of zero 
point. When the sums over all galaxies are taken, these will not average out because the 
direct and inverse magnitudes have different weights. The experience here (Table summed 
to B = 14, Table [7] col. 7 in § 17. 1[ and Table [8] also later in § 17.11) . the mean modulus 
taken over all cluster galaxies is always larger for the direct calculation than for the inverse. 
The modulus difference between the direct and inverse is always smaller when their slope 
difference is smaller, seen from the second half of Table [7J 

In view of these two reasons, which regression should we use to obtain the best distance 
modulus? We have decided not to decide between them but to reduce the data with both, 
avoiding the bias corrections for the direct method by summing to B — 14, a magnitude 
that is symmetrical brighter and fainter than the cross-over point of the direct and inverse 
regressions. 

7.1. The Cepheid and Cluster Zero Point Calibrations Using Various TF 

Slopes 

For each adopted TF slope for the cluster sample we must determine new zero points 
for the Cepheid calibration based on Table [2] of §[2] by forcing the adopted slope on the 
calibrators and calculating the resulting intercept for the calibration equation. (As said 
before, this point is sometimes missed in the literature where a same slope is often used both 
for the direct and inverse method in the Cepheid calibration even as these slopes differ from 
the slopes that apply to the cluster data). 

In the preceding sections we have used a variety of slopes, some arbitrary and others cal- 
culated from the Cepheid data, and some from the composite cluster with various magnitude 
truncations. Table [7| summarizes all of the choices discussed in the previous sections. 
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The equations of the regressions in Table [7] are M = a log W20 + b and B = a log W20 + c 
where the a, b, and c values are listed in columns (1), (2), and (4) of the table. Column (3) 
is the rms of the absolute magnitude zero point calibration, b, in column (2). Column (5) 
is the rms of the apparent magnitude zero point in column (4). Column (6) identifies the 
sample used, ranging from n = 92 for the total sample, n = 84 for the total sample truncated 
at B = 14.0, to n = 78 for a sample still further truncated by eliminating the four outliers 
on the bight side and the two outliers on the faint side near the limit lines in Figures [5] and 
EJ Column (7) shows the mean distance modulus using the cluster data with the listed b 
and c calibrations and the magnitude limits for the cluster data in column (8). The moduli 
in column (7) are the differences between columns (4) and (2). The first five rows are the 
results for the run of assumed slopes discussed in § 16.11 

The strong dependence of the final moduli (summed to the cluster magnitude limit over 
the range of B from about 10 to 15 mag) on the slope of the TF regression is shown by the 
variation in column (7) for the slopes in first five rows of the table. For the more restricted 
range of the slopes shown by the four pairs of values in the last part of the table, the modulus 
values range from 31.30 to 31.60. The values for the highest weight sample (i.e. that with 
the lowest rms values) with n = 78, vary from (m — M) = 31.42 for the inverse solution to 
31.56 for the direct solution. 

It is of interest to set out the data for this highest weight sample in the same detail as 
in Table [5] for the complete sample. Tabled] lists the partial moduli in the various magnitude 
intervals and magnitude limits for the restricted sample format for the direct slope of —6.346 
and the inverse slope of —7.519. Table [HJA. shows the results in magnitude intervals; Table [HJ3 
shows the results of summing to various magnitude limits. 

The data are plotted in Figure [TTJ The open circles are from Table |SJ\; the crosses from 
Table 0B. The noise in the correlations is caused by the nonuniform filling of the correlation 
between B — 12 and 14, mentioned earlier. 



7.2. The Adopted Distance Modulus for the Virgo A, Virgo B and Ursa 

Major I and II Clusters 

From Tableland Figure [TT1 we adopt as our final TF modulus for the composite cluster 
the mean of the direct and inverse calculation using the n = 78 restricted sample. The 
values from Table [33 are (m — M) = 31.56 ± 0.05 for the direct calculation and (m — M) = 
31.41 ± 0.05 for the inverse, giving (m — M) = 31.49 ± 0.04 as the weighted mean. 



Following |FTS 98l (§ 8 there) we apply a systematic correction of 0.07 mag (reducing 
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the modulus) to account for the fact that the cluster members at a given LW are redder 
on average (and also are hydrogen deficient) than for the calibrators. Hence, our final TF 
modulus of the Virgo A cluster is (m — M) = 31.42 ± 0.04 (internal). 

Applying the magnitude offset from Virgo A to Virgo B of 0.38 ± 0.16 mag (Virgo B is 
more distant than Virgo A), of 0.16 ± 0.13 mag for UMa I made closer, and 0.16 ± 0.17 mag 
for UMa II made further from § 14. II gives the modulus values of these clusters as (m — M) = 
31.80 ± 0.16 for Virgo B, 31.26 ± 0.13 for UMa I, and 31.58 ± 0.17 for UMa II. 



7.3. Systematic Errors 

The formal (internal) error of ±0.04 for Virgo A is, of course, too small because of 
systematic errors in (1 ) the Cepheid P-L relation of the order of 0.1 mag in its zero point 



(ITammann et al.ll2007l ). and the error in the zero point of equations §5§ and ([6]) here, (2) the 
uncertainty of ~0.1 mag in the zero point of the Cepheid TF regression in Figure [S] (cf. the 
Cepheid rms value in Table [7]), and (3) the statistical noise due to the nonuniform filling of 
the TF correlation in Figures [5] and [6] as discussed above. It is impossible to put an exact 
value on the systematic accuracy of (m — M) due to these causes, but a reasonable estimate 
is that it is no smaller or larger than ±0.2 mag. Hence, we adopt (m — M) = 31.42 ±0.2 mag 
as our final modulus of Virgo A using the Tully-Fisher method. 

If, as is sometimes done, we combine Virgo A and B, treating them as a single cluster 
with (m — M) = 31.42±0.2 (assumed) from 43 galaxies in Virgo A and (m — M) = 31.80±0.2 
(assumed) from Virgo B from 16 galaxies, the weighted mean of A and B is (m — M) = 
31.52 ± 0.1 4 (extern al assumed), which, however, we do not use. For comparison, the value 
derived by FTS 98l (their eq. 11) where clusters A and B are combined, is 31.58 ± 0.24 
(external) . 



7.4. The Hubble Constant from the Adopted Virgo A Distance 

From the discussion by FTS 98l (Fig. 14 and § 11.2), the expansion velocity of the Virgo 
core, freed from all streaming motions and ti ed to the external cosmic frame b eyond the 
loc al bubble whose edg e is near 4000 km s" 1 flJerien fc TammanrJ Il993l : IFST94 Figs. 17- 
19; iMasters et al.ll2006l ). is 1175kms _1 . This is derived f rom a high weight s o lution of the 
relative Hubble diagram (velocity vs. distance ratios) by iJerjen &: Tammannl (119931 ) where 
the Virgo core is tied to 17 more remote clusters augmented by 24 clusters from Giovanelli 
(private communication in 1996) who also used distance ratios to Virgo. 
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Taking Virgo A to be the core relative to which the distance ratios to the 41 galaxies 
pertains, and using (m — M) = 31.42 ± 0.2 (D = 19.2 Mpc, with a range from 17.5 to 
21.1 Mpc) combined with 1175 km s -1 for the free expansion velocity, the resulting Hubble 
constant is H = 61km s _1 Mpc^ 1 with a range from 53 to 70, where we have also put an 
error of iSOkms^ 1 on the Virgo A expansion velocity. 

If, for some reason, one does not wish to rely on the distance ratios to Virgo of 
Jerjen fe Tammannl (119931 ). as augmented by Giovanelli's private communication in 1996, 
another route to obtain the free expansion velocity of the Virgo core is via a velocity per- 
turbation model for the Virgo velocity relative to the re mote frame. The ob served mean 
heliocentric velocity of the Virgo core is 1050 ±35 km s" 1 fainggeli et alJll993h. This trans- 
form s to 932 km s -1 relative to the ce ntroid of the Local Gro up ( Yahil. Tammann. fe Sandage 
19771 ). Using i>(infail) = 220kms _1 ( Kraan-Kortewegill986l ) for the Virgo infall velocity gives 



irgo 



1152 ± 35kms 1 which is the same as 1175kms 1 to within statistics. 



The value of Hq = 61 obtained here agrees well with the value determined via the super- 
novae rout e in 25 separate stud ies, most of which give H between 50 and 70 (summarized in 
Table 6 of ISandage et al.ll2006l ) by a plethora of authors. However, we consider the method 
here to be of considerably lower weight than that using supernovae because of the evident 
problems raised by the population incompleteness bias when using clusters and the problem 
of tying the Virgo A cluster accurately into the remote cosmic kinematic frame. 
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Fig. 1. — Model of a toy cluster (the insert) with a dispersionless TF correlation of line width 
(corrected for inclination) with apparent magnitude that has a slope of dB^/d log w 2 o = —7. 
A calibration of absolute magnitude with LW in the main body of the diagram is shown with 
variations in the TF slope ranging from —9 to —5. 
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Fig. 2. — (top): The variation of the derived distance modulus at a given sampling magnitude 
into the cluster luminosity function for five assumed TF slope values from data in the first 
half of Table [H All modulus values cross at (m — M) = 31.5 by construction, (bottom): 
The mean modulus summed to the listed cut-off magnitude for the five slope values from 
the second half of Table [U The correct modulus of (m — M) = 31.5 is obtained only when 
B^, 1 stops at 13.0 for all slope values that differ from the true slope of —7. This magnitude 
is symmetrical about the cross-over magnitude of 11.5 (by construction). Sampling fainter 
than the symmetrical magnitude of 13.0 will give incorrect distance moduli, except for the 
case of the correct slope of —7. 
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Fig. 3. — Model of the TF correlation in a cluster with intrinsic dispersion and/or an 
appreciable back-to-front ratio. The bias effect of an incomplete sampling of the cluster 
luminosity function is shown, given, for illustration, an observer's magnitude cut-off limit at 
B^ 1 = 11.5. The right hand vertical LW interval is unbiased because none of the data in 
it are denied entrance into the sample fainter; all are brighter than the cut-off magnitude. 
However, the line width intervals to the left are progressively biased at this magnitude cut- 
off. The incompleteness bias at any given cut-off magnitude is found by summing over all 
line widths brighter than the cut-off value. 
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Fig. 4. — The data for a cluster are binned in apparent magnitude intervals. Four magnitude 
intervals are shown, each enclosed by envelope lines that show the intrinsic dispersion (and/or 
back-to-front effect) in LW at given apparent magnitudes. Slopes for both the direct and 
inverse TF least squares ridge lines for the calibrators are shown schematically. 
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Fig. 5. — The correlation of line width, corrected for inclination, with apparent magnitude 
for the 92 galaxies in the composite cluster composed of galaxies from the Virgo A and B 
subclusters and the Ursa Major cloud as reduced to the distance of the Virgo A subcluster. 
The ridge lines for the direct and inverse least squares regressions are shown whose slopes 
are in parentheses. The dashed envelope lines are pu t one m agnitude brighter and fainter 
than the direct ridge line. The magnitudes are on the IFTS 98l system (their Table 3) which 
is within a few hundredths of a magnitude of the RC3 B^ 1 fully corrected system. The line 
widths are on the LEDA system which is 0.033 ± 0.002 dex smaller than the RC3 system. 
The equations of the ridge lines are equations jT]) and (|2J) of the text. 



-32 - 




2.0 2.2 2.4 2.6 2.8 3.0 
log (W/sin i) (LEDA) 



Fig. 6. — Same as Figure [5] but for a restricted sample of 84 galaxies in the artificial cluster, 
cut from the total sample by B^ 1 < 14.0. The equations of the direct and inverse ridge lines 
are equations (j3J) and (j3J of the text. 
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Fig. 7. — (Top) Change of the Tully-Fisher slope for partial samples as depth of penetration 
into the cluster luminosity function increases using data for the total (n = 92) sample. 
The DIRECT (dots) and INVERSE (open circles) least squares solutions are shown as the 
number of galaxies in the solutions increases from 28 to 92 and the magnitude cut-offs range 
from B^ 1 = 11.0 to 15.0. (Bottom) Same for the restricted sample of 78 galaxies. The data 
are from Table O 
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Fig. 8. — The Cepheid calibrator data from Table [3] shown with the least squares regressions 
for both the direct and inverse formulations whose slopes are in parentheses. These differ 
from the slopes from the data for the composite cluster in Figures and El both of which 
have higher weight than the Cepheid sample alone. 
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Fig. 9a. — The variation with apparent magnitude of the calculated distance moduli for each 
galaxy in the artificial cluster when binned in logu^o intervals in steps of 0.1 dex using the 
DIRECT slope of —6.622 in Table H]A. for the complete sample of n = 92. The mean modulus 
marked by the arrow in each panel is that listed at the bottom of Table 0]A. A light vertical 
line to guide the eye is put in each panel at B — 12.0. 
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Fig. 9b. — Same as Figure [9^ using the INVERSE slope of —8.032 and the individual data 
for the composite cluster. 



-37- 



cn 
_i 

Z> 
Q 

o 



32.2 
32.0 
31.8 
31.6 
31.4 
31.2 

31.8 
31.6 
31.4 
31.2 



DIRECT 

-6.622 
n = 92 



n = 92 




o- between B°-'_ 
*- up to B°j' 



+ 



+ 



-i 1— 



INVERSE 

-8.032 



n = 92 



10 11 12 13 14 15 



Fig. 10. — (Top). The average moduli for the total sample (n = 92) listed in Table [5] in 
each magnitude interval using the direct slope of —6.622 and summed over all line widths 
in that magnitude interval (open circles), and then summed over all magnitude intervals 
(crosses) up to the listed Bj? magnitude cut-off limit (crosses). (Bottom) Same for the 
inverse formulation with a slope of —8.032 for the total sample. The noise in both the open 
circle and the crosses data between B = 12 and 13 is caused by nonuniform filling of the TF 
distribution between the envelope lines of Figures [5] and El The data are from Table [5j 
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Fig. 11. — Similar to Figure [TU1 but for the restricted sample with n = 78. The Teerikorpi 
cluster population incompleteness bias is present using the direct TF slope but is absent 
using the inverse slope. The open circles are the modulus values from the data summed 
over all LW and averaged within the individual magnitude intervals (Table [SJA.). The crosses 
are the modulus values for the data summed aver all LW and then summed to the various 
magnitude cut-off limits (Table HP). 
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Table 1. Derived Distance Moduli of the Toy Cluster of Figured] as a Function of B^ 1 
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Table 2. Least Squares TF Regressions to Various Magnitude Cut-offs for the Total 
(n = 92) and the Restricted (n = 78) Sample of the Composite Cluster 
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DIRECT INV DIRECT INV 
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Table 3. The Cepheid Calibrators 



Name 


Type 
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log W 20 
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Sbl-II 


24.54 


3.36 


-21.18 


2.728 


NGC300 


ScII.8 


26.48 


8.49 


-17.99 


2.296 


NGC598 


Sc(s)II-III 


24.64 


5.73 


-18.91 


2.357 


NGC925 


SBc(s)II-III 


29.84 


10.04 


-19.80 


2.382 


NGC1365 


SBbc(s)I 


31.46 


9.94 


-21.52 


2.648 


NGC 1425 


Sb(r)II 


31.96 


10.83 


-21.13 


2.582 


NGC2090 


Sc(s)II 


30.48 
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NGC 3319 
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Table 4. (m — M) Averaged within Various Discrete Intervals of Line Width and Listed in 

Magnitude Intervals for the Composite Cluster 
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32 50 (11 


31.79(2) 


31.61(1) 


14.0-15.0 












32.82(1) 


32.12(7) 


Summed 
















vertical 


31.60 


31.72 


31.43 


31.58 


31.65 


31.42 


31.77 


n 


5 


10 


16 


17 


20 


14 


10 


rms 


0.145 


0.563 


0.483 


0.372 


0.556 


0.573 


0.667 






(B) INVERSE SLOPE = - 


-8.032 






9.8-10.5 


31.85 (5) 


31.29 (3) 


30.56(1) 










10.5-11.0 




31.79 (3) 


30.93 (4) 










11.0-11.5 




31.88 (2) 


31.32(7) 


30.98 (2) 


30.29(1) 






11.5-12.0 




32.68 (2) 


31.86 (4) 


31.36 (6) 


30.51(2) 






12.0-12.5 








31.60 (6) 


30.13(1) 






12.5-13.0 








31.99 (2) 


31.47(5) 


30.61(6) 


29.69(1) 


13.0-13.5 










31.84(5) 


31.26 (4) 


30.33(1) 


13.5-14.0 










32.19(1) 


31.34(2) 


31.05(1) 


14.0-15.0 












32.41(1) 


31.54(7) 


Summed 
















vertical 


31.85 


31.84 


31.41 


31.43 


31.27 


30.99 


31.18 


n 


5 


10 


16 


17 


20 


14 


10 


rms 


0.179 


0.552 


0.481 


0.389 


0.561 


0.580 


0.673 
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Table 5. (m — M) as Summed over all Line Widths for Various Magnitude Intervals and 

Cut-off Magnitudes for the Composite Cluster 



between 




summed over LW 


up to 




summed 


over LW 


nOi 

JDrp 


n 


(m - M) 


rms 


R 0,« 

JDrp 


n 


(m-M) 


rms 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 






(A) DIRECT SLOPE 


= -6.622 








9.8-10.5 


9 


31.32 


0.398 


10.5 


9 


31.32 


0.398 


10.5-11.0 


7 


31.25 


0.470 


11.0 


16 


31.29 


0.417 


11.0-11.5 


12 


31.45 


0.410 


11.5 


28 


31.36 


0.414 


11.5-12.0 


14 


31.60 


0.570 


12.0 


42 


31.44 


0.479 


12.0-12.5 


14 


31.44 


0.427 


12.5 


56 


31.44 


0.463 


12.5-13.0 


14 


31.40 


0.550 


13.0 


70 


31.43 


0.478 


13.0-13.5 


10 


31.82 


0.402 


13.5 


80 


31.48 


0.484 


13.5-14.0 


4 


31.92 


0.413 


14.0 


84 


31.50 


0.488 


14.0-15.0 


8 


32.21 


0.286 


15.0 


92 


31.56 


0.514 






(B) INVERSE SLOPE 


= -8.032 








9.8-10.5 


9 


31.52 


0.475 


10.5 


9 


31.52 


0.475 


10.5-11.0 


7 


31.30 


0.550 


11.0 


16 


31.42 


0.504 


11.0-11.5 


12 


31.40 


0.510 


11.5 


28 


31.42 


0.495 


11.5-12.0 


14 


31.50 


0.693 


12.0 


42 


31.44 


0.561 


12.0-12.5 


14 


31.20 


0.527 


12.5 


56 


31.38 


0.558 


12.5-13.0 


14 


31.05 


0.676 


13.0 


70 


31.32 


0.594 


13.0-13.5 


10 


31.44 


0.495 


13.5 


80 


31.33 


0.581 


13.5-14.0 


4 


31.48 


0.300 


14.0 


84 


31.34 


0.576 


14.0-15.0 


8 


31.65 


0.354 


15.0 


92 


31.36 


0.566 
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Table 6. (m — M) as Summed to Various Magnitude Cut-offs over all LW for Different 

Tully-Fisher Slopes 



Summed to SLOPES 

-9 -7 -6 -5 

By 1 n (m — M) rms (m — M) rms (m — M) rms (m — M) rms 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



10.5 


9 


31.66 


.524 


31.37 


.418 


31.23 


.365 


31.09 


.317 


11.0 


16 


31.51 


.567 


31.33 


.439 


31.23 


.383 


31.14 


.339 


11.5 


28 


31.45 


.563 


31.37 


.433 


31.33 


.388 


31.29 


.366 


12.0 


42 


31.44 


.632 


31.44 


.498 


31.44 


.455 


31.44 


.436 


12.5 


56 


31.34 


.643 


31.43 


.485 


31.46 


.437 


31.50 


.423 


13.0 


70 


31.23 


.697 


31.40 


.504 


31.48 


.447 


31.56 


.434 


13.5 


80 


31.23 


.678 


31.44 


.504 


31.55 


.467 


31.65 


.497 


14.0 


84 


31.22 


.671 


31.46 


.505 


31.57 


.477 


31.69 


.497 


15.0 


92 


31.23 


.651 


31.51 


.518 


31.65 


.525 


31.78 


.576 
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Table 7. Regressions for the Cepheid Calibration and for the Composite Cluster for 
Various Tully-Fisher Slopes: the Sensitivity of the Resulting Distance Moduli to Slope 



blopc 


Ccphcid calib. 5 ) 


Cluster 6 ) 


n 


(to — 


M) 


R 0,i 


a 


b 


rms 


c 


rms 




to sample end 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


-9 


+2.716 


0.624 


33.946 


0.651 


92 


31 


23 


15.0 


-8 


+0.197 


0.541 


31.525 


0.564 


92 


31 


33 


15.0 


-7 


-2.402 


0.511 


29.108 


0.518 


92 


31 


51 


15.0 


-6 


-4.961 


0.499 


26.689 


0.525 


92 


31 


65 


15.0 


-5 


-7.519 


0.522 


24.261 


0.576 


92 


31 


78 


15.0 


-6.161 1 ) 


-4.549 


0.500 


27.078 


0.520 


92 


31 


63 


15.0 


-8.475 1 ) 


+1.364 


0.580 


32.664 


0.598 


92 


31 


30 


15.0 


-6.622 2 ) 


-3.369 


0.502 


28.191 


0.514 


92 


31 


56 


15.0 


-8.032 2 ) 


+0.239 


0.557 


31.609 


0.566 


92 


31 


37 


15.0 


-5.936 3) 


-5.124 


0.500 


26.456 


0.477 


84 


31 


58 


14.0 


-7.547 3 ) 


-1.003 


0.533 


30.387 


0.538 


84 


31 


39 


14.0 


-6.346 4 ) 


-4.077 


0.499 


27.483 


0.415 


78 


31 


56 


14.0 


-7.519 4 ) 


-1.077 


0.530 


30.343 


0.451 


78 


31 


42 


14.0 



1 - ) Cepheid least squares regressions, n = 25 

2 ) Direct and inverse slopes for total cluster sample (Fig. [5]) 

3 - ) Direct and inverse slopes for cluster sample cut at B^ 1 — 14.0 
(Fig.® 

4 ) Direct and inverse slopes for final sample (with a B = 14.0 cut plus 
six others eliminated) 

5 ' Cepheid calibration as Mbit) = fl logW20 + b 

6 ) Cluster regression as B^ 1 = a log W20 + c 
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Table 8. Modulus of the Composite Cluster for Various Magnitude Intervals and Limits 

from the High Weight Restricted Sample 







Slopc= - 


-6.346 


Slopc= - 


-7.519 


Interval 




Direct (78) 


Inverse (78) 


R 0,i 

t>rp 


n 


(to - M) 


rms 


(m - M) 


rms 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(A) IN MAGNITUDE INTERVALS 


9.8-10.5 


8 


31.78 


.270 


31.56 


.311 


10.5-11.0 


7 


31.27 


.505 


31.28 


.521 


11.0-11.5 


11 


31.53 


.299 


31.52 


.355 


11.5-12.0 


12 


31.47 


.407 


31.35 


.482 


12.0-12.5 


13 


31.54 


.339 


31.36 


.412 


12.5-13.0 


13 


31.55 


.446 


31.28 


.535 


13.0-13.5 


10 


31.90 


.383 


31.58 


.461 


13.5-14.0 


4 


32.01 


.396 


31.64 


.480 


(B) TO MAGNITUDE CUT-OFF LIMITS 


10.5 


8 


31.38 


.270 


31.56 


.311 


11.0 


15 


31.33 


.386 


31.36 


.442 


11.5 


26 


31.42 


.360 


31.43 


.407 


12.0 


38 


31.43 


.371 


31.41 


.427 


12.5 


51 


31.46 


.363 


31.39 


.420 


13.0 


64 


31.48 


.379 


31.37 


.444 


13.5 


74 


31.54 


.404 


31.40 


.449 


14.0 


78 


31.56 


.415 


31.41 


.451 



